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for study on finite- (T,/_f) QCD since its discrete symmetries are appropriate for the 
case. To show its applicability, we study strong-coupling lattice QCD with temper- 
ature and chemical potential. We derive the effective potential of the scalar meson 
field and obtain a critical line of the chiral phase transition, which is qualitatively 
consistent with the phenomenologically expected phase diagram. We also discuss 
that 0{l/a) renormalization of imaginary chemical potential can be controlled by 
adjusting a parameter of a dimension-3 counterterm. 
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I. INTRODUCTION 

Understanding of QCD(Quantum Chromo Dynamics) under extreme conditions with 
temperature and density is one of keys to elucidating the history of the universe. In partic- 
ular, the QCD phase diagram has been attracting a great deal of attention (See for example 
[l|). Lattice QCD is the most powerful tool to investigate such non-perturbative aspects of 
QCD. Indeed, the lattice QCD simulation has been applied to the finite-temperature QCD 
with zero density, and has produced lots of works to investigate the critical or crossover be- 
havior due to confinement and deconfinement transition (See for example However, the 
Monte-Carlo simulation cannot be easily applied to QCD with chemical potential because 
of the notorious sign problem (See references in jsj). There have been developed several 
prescriptions to bypass this problem, including the imaginary chemical potential method, 
ihe Taylor expansion, the Fugacity expansion and the histogram method (See references in 
^). Apart from the numerical simulation, the analytical lattice study have been also devel- 
oped. One of the classical and reliable methods is the strong-coupling expansion 0-0]. This 
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lod has been applied to the QCD phase diagram and has produced successful results 



18| . In these works, (unrooted) staggered fermions 19l-l22| have been used, thus the 
corresponding continuum theory is 4-fiavor QCD although the physical two or three-flavor 
QCD are desirable. 

In this paper we propose a new framework of investigating the 2-fiavor finite- (T, n) QCD 



phase diagram by using the Karsten- Wilczek (KW) lattice fermion discretization [23| . This 
lattice formulation lifts degeneracy of 16 species by introducing a species-dependent (imag- 
inary) chemical potential term, instead of introducing a species-dependent mass term in 
Wilson fermion. The most notable point is that it can reduce the number of species to 
2 with keeping U{1) chiral symmetry. It is sometimes called "minimal-doubling fermions" 



23N25| or "flavored-chemical-potential(FCP) fermions" [26|. The phase structure in the pa- 



rameter plane for them has been recently studied in 26|]. In the present work we show that 



the KW discretization suits study on the 2-flavor finite-(T, /i) QCD phase diagram since it 



has the same discrete symmetries 27H29| as the finite-density lattice QCD system. With 
progress on the sign problem, this formulation can be a powerful tool for the in-medium 
lattice QCD. To show the usefulness of the KW fermion, we study strong-coupling lattice 
QCD with temperature and density. We derive the mesonic effective potential as a function 
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of (T,/i) and elucidate a critical line of the chiral phase transition. The result is qualita- 
tively consistent with predictions from the phenomenological models. Toward a practical 
application to lattice simulations, we also argue that the 0{l/a) renormalization of imag- 
inary chemical potential can be controlled by adjusting the relevant parameter /i3 of the 
dimension-3 operator. 

In Sec. [TTl we investigate the KW fermion and argue that it is a useful formulation for 
two-flavor (T,yu) lattice QCD. In Sec. IIIII we study the strong-coupling lattice QCD and 
derive the QCD phase diagram. Section IIVI is devoted to a summary and discussion. 



II. FLAVORED CHEMICAL POTENTIAL 



The Karsten-Wilczek (KW) fermion discretization decouples 14 among the 16 species in 
the naive fermion by introducing a species-dependent ima gina ry chemical potential term 
without losing all the chiral symmetry and ultra-locality 23[. This is a special case of 



"Flavored-chemical-potential (FCP) fermions'' 



which includes a real-potential type and 



an imaginary-potential type. The KW fermion is an imaginary-type FCP fermion and can 



-doubling" 
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23|-|25|. 



2811. we need to 



describe two flavors with a proper parameter value, as is called "minima 
Since the flavored chemical potential term breaks discrete symmetries 
flne-tune the three parameters for one dimension-3 (tpi'-f^il)) and two dimension-4 {ipd^^p, 
Fj^^Fjij counterterms in order to take a correct Lorentz-symmetric continuum limit for the 
zero-(T,/i) QCD simulations 29|. However, as we will see later, the discrete symmetries of 
KW fermions suit the flnite-temperature and -density system, and the severe flne-tuning to 
restore Lorentz symmetry would not be required for this case: What we need to care about 
in this case is the dimension-3 operator which corresponds to an 0(l/a) chemical potential 
term. This fact inspires us to apply the KW fermion to the in-medium QCD. 



A. Symmetries 

We study the Karsten-Wilczek fermion and its symmetries in comparison with other 
lattice fermions with chemical potential. The 4-d KW action with f/(l) chiral symmetry 
and ultra-locality is obtained by introducing a Wilson-like term proportional to ^74 into the 
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naive fermion action as 



KW 



E 



+i/i3?/;^74V'x 



(1) 



where the second hne with a parameter r is a flavored-chemical-potential term, which works 
to hft the degeneracy of species. The third hne with a relevant parameter p,^ is a dimension-3 
counterterm, which corresponds to an 0(1 /a) chemical potential term 4^. If one drops i 
in front of the second and third lines, this becomes a real-type Karsten-Wilczek fermion 
without 75 hermiticity while we in this paper focus on the imaginary-type KW fermion 
basically. For the free theory, the associated Dirac operator in momentum space is given by 

4 3 

aDKw (p) = 7m sin ap/. + ^74 (/^s 3r - r ^ cos apj). (2) 
M=i j=i 

For r = 1 and fi^ = 0, it has only two zeros at p = (0, 0, 0, 0), (0, 0, 0, vr/a). The rest 
14 species have 0(1/ a) imaginary chemical potential due to the flavored-chemical-potential 
term. More precisely, among the original 16 species, two have zero imaginary chemical 
potential, six have 2/a, six have A/a and two have 6/a. In the naive continuum limit, the 
14 species are decoupled with infinite imaginary chemical potential and there remains only 
two flavors as shown in Fig. [D j45|. The flavored chemical potential term breaks discrete 
symmetries. The residual symmetries are cubic symmetry, corresponding to permutation of 
spatial three axes, CT and P 27|]. We list the symmetries of importance as following; 



(1) U{1) chiral symmetry (75 

(2) P 

(3) CT 

(4) Cubic symmetry. 



T3 



30 



Since the cubic symmetry is likely to be enhanced to the 3d rotation symmetry in the 
continuum, we expect that these symmetries become those of the finite-density QCD in the 
continuum limit. To convince ourselves, let us look into symmetries of the naive lattice 
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FIG. 1: Species-splitting in Wilson and Karsten-Wilczek fermions. Circled numbers stand for the 
number of massless flavors on each point. 

fermion with chemical potential, which is introduced as a 4th-direction abelian gauge field. 

The massless naive action with chemical potential is given by 

■ 3 

^ ] '4'x'Jj iUx,x+j4'x+j ~ Ux,x-j'4'x-j) 



(3) 



The chemical potential breaks the hypercubic symmetry into the spatial cubic symmetry. 
It also breaks C, P and T symmetries into CT and P. These discrete symmetries are the 
same as those of Karsten-Wilczek fermion as shown above. From the viewpoint of the 
universality class, these two theories belong to the same class. We therefore consider that 
the KW fermion suits in-medium lattice QCD much better than zero-(T,/i) lattice QCD. 

However, we have to care about the way of introducing chemical potential. In KW 
fermions the flavored imaginary chemical potential is introduced naively as ipi'y^ip while the 
chemical potential is usually introduced as a 4th-direction abelian gauge fleld on the lattice as 
shown in ([3]). As is well-known 35|, the naive introduction of chemical potential violates the 
abelian gauge invariance and requires a counterterm to make the energy density and other 
thermodynamical quantities flnite. We can discuss the necessity of a counterterm also from 
the viewpoint of the additive renormalization: We remind ourselves that the Wilson fermion 
breaks chiral symmetry by 0(1/ a) flavored-mass terms, which leads to the additive mass 
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renormalization and the necessity of the mass parameter tuning in the interacting theory. 
Now the KW fermion breaks discrete symmetries into those of finite-density systems by 
0(l/a) flavored-chemical-potential terms. These facts indicate that we will here encounter 
large chemical potential renormalization instead of the additive mass renormalization in 
Wilson fermion. This is why we need to introduce the dimension-3 counterterm as fisipijii/j 
in the action even for the application to finite-density QCD. The zero-chemical-potential 
two flavors for a free case in ([2]) suffer 0(1/ a) imaginary chemical potential renormalization 
in the interacting theory. We thus need to tune fi^ to control it even if we apply it to the 
finite-density QCD. 



B. Additive chemical potential renormalization 



As the additive mass renormalization in Wilson fermion is manifested in the phase dia- 
gram in a {m-g^) plane 36l-l40l] . the chemical potential renormalization can be manifested 
in a phase diagram in a (/is-f?^) plane. The chiral phase structure in the parameter space of 
lattice QCD with the KW fermion is studied in [26| by using strong-coupling lattice QCD 
and the Gross-Neveu model: Fig. [2] is the conjectured chiral phase diagram with the num- 
ber of physical flavors in the {^i^-g"^) plane for r = 1. There are roughly two phases with 
and without chiral condensate, or equivalently with and without SSB of chiral symmetry. 
We name them as "physical" and "unphysical" phases since the physical QCD has SSB of 



chiral symmetry at least. As shown in 
boundaries between the two phases are given by 

6r2 + 2 



261], in the strong-coupling and large limits, the 



/is 



— 3r, 



(4) 



which gives the physical range — ^32/7 — 3 < /ia < + ^^32/7 — 3 for r = 1 . We thus have 
the two chiral boundaries in this limit as shown in Fig. [2l 

In the weak-coupling limit {g"^ = 0) we analytically know the number of flavors. In 
([2]) with r = 1, the number of flavors changes with fi^ being varied 2^'■ There are four 
sectors with two, six, six and two flavors as shown in Fig. [2l (On boundaries between the 
sectors fi^ = 1,-1, —3, —5, —7, fermions have unusual dispersions as D{p) ~ p + p^, which 
cannot be fixed even by tuning parameters. We therefore avoid these points.) There are 
no fermion flavors outside these four sectors, but only unphysical fermions with 0(1/ a) 
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chemical potential. 

As seen from Fig. [2l the boundaries between physical (a) ^ and unphysical (cr) = 
phases start from boundaries between the two-flavor and no-flavor sectors in the weak- 
coupling limit. It is reasonable since SSB of chiral symmetry can take place only in theories 
with fermions. We especially call the two-flavor range "minimal-doubling range" . For r = 1 
the minimal-doubling range is given hj —1 < fi^ < 1 and —7 < fi-s < —5. Toward the strong 
coupling limit, these ranges are expected to change with g"^ as Fig. [2j We note that the 
minimal-doubling range and the 6-flavor range become less distinguishable with the gauge 
coupling being larger. 

From the viewpoints of practical application, the relevant parameter fi^ has to be tuned 
to cancel the 0(1/ a) imaginary chemical potential renormalization for the two flavors as 
discussed in the previous subsection. In the weak-coupling limit, it is obvious that the two 
flavors feel no chemical potential as long as we set fi^ in the minimal-doubling range as 
— 1 < /i3 < 1 or —7 < /i3 < —5 for r = 1. We consider that it carries over in the interacting 
theory: For a given gauge coupling, in order to cancel the 0(1/ a) effective chemical potential 
for the two flavors, we have to set a value of fi^ in the minimal-doubling range. As conjectured 
in Fig. 121 the minimal-doubling range for the middle gauge coupling is likely to have some 
width, which means that we do not need to fine-tune /is, just set it in the range. After we 
tune /i3 as such, there will remain only 0(1) renormalization of imaginary chemical potential 
which has a physical scale. 

From these arguments, it becomes clear that KW fermion can be applied to study on the 
finite- (T,/i) QCD as long as /is is chosen to be a proper value. We here write the KW action 
with a usual chemical potential parameter /i as 



where /i is a chemical potential parameter, which can be real /iRe, imaginary ifiim or complex 
/iRe + ^Aiim- As we discussed, we consider that we can keep the additive 0(1/ a) renormaliza- 
tion of imaginary chemical potential under control by setting /is within the range. We also 
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+ V^x74 C^lpx - U^,x+j1px+j - U^,x-ji^x-j) + ifl'si^xlii^x 

i=i 



(5) 
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FIG. 2: Conjecture on the f^s-g"^ chiral phase structure for the Karsten-Wilczek fermion with r=l. 
The width of the minimal-doubling range determines how hard it is to tune fj,^ 26l |. 



need to adjust /iim to control the 0(1) imaginary chemical potential. We expect we can be 
informed of the size of the effective imaginary chemical potential by condensate in small 
T and /iRc region as we will discuss later. 

In the end of this section, we give a comment on practical applications of KW fermions to 
numerical simulations. Since the present lattice simulation can be applied only to systems 
with small /irc or with /iim, the most feasible application of the KW fermion for now could be 
the imaginary-chemical-potential lattice QCD. In such a case we can in principle describe two 
flavors with arbitrary fii^ by keeping fi-^ within the minimal-doubling range and adjusting 
/iim properly. 



III. STRONG-COUPLING LATTICE QCD 



To show that the Karsten-Wilczek formulation works in the study of finite-temperature 
and finite-density lattice QCD, we study the QCD phase diagram in the framework of strong- 
coupling lattice QCD with this formulation. The strong-coupling lattice QCD study with 
minimal-doubling fermions has been first performed in 4]J], where spontaneous chiral sym- 
metry breaking due to chiral condensate is observed. We extend this to finite-temperature 
and finite-density cases by following a method in 10l4l2l|. and elucidate the dependence of 



chiral condensate on temperature and chemical potential. We will find a phase structure 
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consistent with the phenomenological models. 

Before starting analysis, we give some comments on our analysis: Firstly, in the strong- 
coupling limit, the notion of "species" gets ambiguous as shown in Fig. [2j We cannot 
distinguish 2-flavor and 6-flavor ranges in the {fis,g^) plane, but we can just distinguish 
physical {{a) ^ 0) and unphysical {{a) = 0) regions in the strong-coupling. In this section 
we will just choose a value of fi^ within the physical range, and obtain the QCD phase 
diagram. Our purpose here is just to show that the KW formulation works to study finite- 
(T, /i) QCD, thus the analysis with this modest condition is sufficient. Secondly, as we have 
discussed, we consider the feasible numerical application of the KW fermion would be the 
imaginary-/! lattice QCD for now. However, we in this section introduce real-// and study 
the finite- (T, /irc) phase structure since the strong-coupling lattice QCD is free from the 
sign problem and the results can be compared to the phenomenological predictions. 



We now derive the effective potential of the scalar meson field a from ([5]) in the strong- 
coupling limit [g^ — )■ oo). We here consider general color number N^. for the SU{Nc) gauge 
group and general space-time dimensions a.sd = D-\-l. We first perform the 1-link integral for 
the gauge field in the dimensional spatial part, and introduce auxiliary fields to eliminate 
the 4-point interactions as 



A. Effective potential 




X j = l 



D 




with 
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where we introduce the mesonic field as 



(8) 



We note that two auxihary fields a and 714 are required to get rid of four-fermi interactions 
in this case. 714 condensate is related to density, and we will discuss it later. We also note 
that we dropped the next-leading order of 0(1/ VD) expansions in ([S]), which corresponds 
to a large D limit. We now have an intermediate form of the effective action as 



5'eff - Yl 



— r27r474 ) M.\ 



X) 



(9) 



We here consider real chemical potential as /i = /iRe- We perform Fourier transformation of 
the temporal direction by introducing Matsubara modes as, 



1 



2tx 



n 



' n=l * n=l ^ 

We here take the Polyakov gauge. The link variable in the temporal direction is given by, 

/ „i(Pl{x)/Nr \ 



(10) 



UJx) 



J(f>2(x)/Nr 



\ 



J]0a(5') = O, 



(11) 



a=l 



J 



with (pa defined as components of gauge fields. It enables us to calculate fermionic determi- 
nant analytically as. 



detD = nnn^^* 

X a=ln=l 



D 



m 



D 



+ — Vl + r2(T 1 + 274 sin A;^") + fi^ + Dr - —Vl - t^tta 



X a=l n=l 
X a=ln=l 

Nc 

Yl A'^^^^^ Yl (2 cosh N^E + 2 cos (0„ - iN^fi))' 



(12) 



a=l 
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where we define 



= 1 + + Dr - y Vl - r27r4^ 



D 



B = m -\- — Vi + r^a, 



with kr^^ = kn + 4>a/Nr — ifi, and k^' is determined by the relation Asink^' = sin kr + 



E = arcsinh ( — 
A 



log 



B 



- + W1+ - 



A 



B 



A 



(a) 



(13) 



(14) 



/i3 + Dr — -jVl — T^7r4. By integrating the temporal gauge field (pa we derive 



det ((5n+j-j)i<jj<^^ 



„ AT, 

/ J]^ JJ (2 cosh NrE + 2 cos (0^ - «iVr/i))^ = n 

X a=l X 

(15) 

Q^= [ ^(2coshiV^E + 2cos^)^e-*"'^, 9 = cj) - iN^^i. (16) 
For A'^c = 3 these Q„ are explicitly given as 

Qo = 2(8 cosh^ NrE + 24 cosh^ NrE + 3), Q±i = 8 cosh NrE{A cosh^ iV^E + 3)e^^^'^, 
g±2 = 4(6 cosh^ iV^E + l)e^2^-'^, g±3 = 8 cosh NrE e^^^^^", 



Q±4 = e±^^^^ Q|„ 



l>5 



0. 



(17) 



As a result, the effective potential is given by 

NrD 



4 



((1 + r^)a^ + (1 - r^)7rl) - N, log A 



--log 5^ det (g n+i-j, 

\neZ 



(18) 



Here we redefine the free energy AJ^cS — >■ J^cs to be consistent with the phenomenological 
result as discussed later. We here show only the calculation result of the determinant part 
for A^, = 3, 



5^ det (Q„ 



+i-jJl<i,j<Nc 



E 



E 



T 



T 



E 



E 



E 



1 + 12 cosh^ — + 8 cosh* — 15 - 60 cosh^ — + 160 cosh* — - 32 cosh" — + 64 cosh 



T 



T 



E 



+80 cosh 



T 
T 



T 



E 



E 



E 



64 cosh ^ cosh — - 15 + 40 cosh^ 77; + 96 cosh* — + 320 cosh^ — 



T 



T 



T 



E E E 

1 + 6 cosh^ — + 24 cosh^ — + 80 cosh^ — 
T T T 



3/i 



B 



E 



E 



+80 cosh — — cosh — — 1 + cosh — + 2 cosh 



4^ 
T 
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with 



(20) 



In the case of zero temperature T = 0, we can solve the equihbrium condition analytically. 
For D = 3 {d = 4) with m = and r = 1 the free energy is given by 



J^csi<y, T^A] m, T, fi, /ig) = ^0-^ - ^ log (l + (/is + 3)^) 



max < 3 arcsinh 



3v/2, 



(21) 



2^1 + {^is + S) 

In this case there are two local minima of the free energy as a function of a 
J^cs = -/is - I log (1 + (/is + 3)2) at a = and T^s = - | log (1 + (/is + 3)^) - 



3arcsinh(3 cr/A/2(l + (/is + 3)^)) at cr = ctq. This (Tq satisfies the following gap equation, 



eff 



da 

Therefore we have 



2a' 



<T=ao 



an 



9 



an 



21 + (/is + 3)2 



l + (/^3 + 3)2- 



(22) 



1 + (/i3 + 3)- 
9 



9 



1 



(23) 



(l + (/i3 + 3)T 

Comparing these two local minima, we can show that the global minimum changes from 
(T = ctq to (T = at the critical chemical potential as 

3V2ao 



^criticai(^y = 0) = 3 arcsinh 



2v/l + (/is + 3)^ 



- -<Jn 



(24) 



This phase transition is of 1st order because the order parameter a changes discontinuously 
at this critical chemical potential. We can also evaluate the baryon density Pb = —dJ^cs/ djis 
at T = 0. It turns out to be empty Pb = ^ when Pb < p'b^^^^^^- On the other hand, when 
Pb > p%^^^'^^^, it is saturated as pb = 1. 



B. Phase diagram 

We depict the QCD phase diagram with respect to chiral symmetry from the effective 
potential f llSp . We first concentrate on the case with r = 1 for simplicity. In this case the free 
energy is independent of ir^ although it is an artifact of the strong coupling limit. Thus we 
simply neglect 114 and set ps in the physical parameter range (— -\/32/7 < /is + 3 < +^^32/7) 
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FIG. 3: Phase diagram for the chiral transition with r = 1, fj,^ = —0.9 and m = 0. Green and red 
Unes show 2nd and 1st transition hnes, respectively. The transition order is changed from 2nd to 
1st at the tricritical point (;ug',r*") = (0.402,0.119). 

in the strong-coupling limit in Fig. |2l We here take /is = —0.9. In a massless case m = 0, 
the phase boundary of the 2nd order chiral phase transition is given by the condition, such 
that the coefficient of cr^ in the free energy becomes zero. When the order of the phase 
transition is changed from 2nd to 1st, the coefficient of cr^ as well as should vanish in the 
free energy J'cff(a-). 

Figure [3] shows the phase boundary of the chiral transition with r = 1, = —0.9 and 
m = 0. The order of the phase transition is changed from 2nd to 1st at the tricritical 
point (/i§\T*") = (0.402,0.119). We also depict a condensate and the baryon density 
Pb = —dJ^cff/dpB as functions of ps with several fixed T in Fig. HI We find that there are 
first (T < T*") and second (T > T*''') order phase transitions for cr, followed by the phase 
transition of the density ps- For m 7^ 0, we can easily show that the crossover transition 
instead appears with the 2nd-order critical point. 

These results are qualitatively consistent with those with strong-coupling lattice QCD 
with staggered fermions, while there are some quantitative differences. For example, the 
KW phase diagram is suppressed in T direction compared to that in staggered. We here 
compare the ratio of the transition baryon chemical potential at T = to the critical 
temperature at hb = 0, i?" = fic{T = 0)/Tc{pb = 0). In staggered fermion, this ratio is 
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FIG. 4: Chiral condensate a and the baryon density pB for (left) T = 0.15 and (right) T = 0.10. 
Top and bottom panels show the massless m = and massive m = 0.1 cases. There are 1st and 
2nd phase transitions for a. In the case of m 7^ 0, there appears the crossover behavior instead of 
the 2nd order transition. 



■st 



3 X 0.56/(5/3) 
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12], while - 0.406/0.178 ~ 2.3. In the real world, 



this ratio is larger, ^ Mn/170 MeV ~ 5.5. When the finite coupling and Polyakov 
loop effects are taken into account for staggered fermion, Tc{fiB = 0) decreases, /ic(T = 
0) stays almost constant, then value increases 14l-ll6|. Larger with KW fermion 
in the strong coupling limit may suggest smaller finite coupling corrections in the phase 
boundary. Another interesting point is the location of the tri critical point. In KW fermion, 
the ratio i?Kk,= 0.402/0.119 ~ 3.4, while Rlf = 1.73/0.866 ~ 2.0 for unrooted staggered 



fermion 
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12l |. It would be too brave to discuss this value, but -Rkw consistent with the 
recent Monte-Carlo simulations 42|], which implies that the critical point does not exist in 
the low baryon chemical potential region, Hb/T ^ 3. These observations reveal usefulness 
of KW fermion for research on QCD phase diagram. 

The dependence of a and pb seems to show there are two sequential transitions with 
increasing /i^. AtT = 0.15> T*", a quickly decreases and pB increases at pb — 0.25 
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FIG. 5: Phase diagram for the chiral transition with r = 0.75, ^3 = —0.9 and m = 0. The transition 
order is similarly changed from 2nd to 1st at the tricritical point (/i^^T*''') = (0.450,0.133). 
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FIG. 6: Condensates a and |7r4| for (left) T = 0.15 and (right) T = 0.10 in the case of r = 0.75. 
There are 1st and 2nd phase transitions for a. 

(0.41) for m = (0.1), and at a larger {fJ'B — 0.64), increasing rate of ps as a function 
of fiB becomes higher again. At lower temperature, T = 0.10 < T*"^', partial restoration of 
the chiral symmetry is seen before the first order phase transition. Since we have not taken 
care of the diquark condensate, these transitions are not related to the color superconductor. 



431] , partial chiral restored matter 



Other types of matter, such as quarkyonic matter 
or nuclear matter, may be related to the above sequential change. 

We next consider r 7^ 1. We as an example take r = 0.75 and fis = —0.9, which is 
again within the physical range. (For r = 0.75, the physical parameter range is given by 
—3.84 ^ /is ^ —0.66 from (j4]).) Fig. |5] shows the phase diagram for the chiral transition. 



16 



which is similar to that of r = 1 case. We have a httle higher critical temperature, Tc{fiB = 
0) ~ 0.20, and a little larger transition chemical potential, fic{T = 0) ^ 0.42, compared 
with r = 1 case. The tricritical point, (/i^\T*") ~ (0.450,0.133), also has a larger T and 
fiB- One difference is the existence of the region, where dfis/dT > on the first order 
transition boundary. This is a lattice artifact, and can be removed when we take account 
of finite coupling effects. The Clausius-Clapeyron relation tells us that the balance of the 
pressure in the coexiting two phases, I and II, leads to the slope of the phase boundary, 
d^s/dT = —(su — Si)/(pn — Pi), where Sj and pi represent the entropy and baryon density 



in the phase i [12[. The positive slope of the phase boundary thus implies that the entropy 
density in quark matter (higher fiB phase) is smaller than that in hadronic matter (lower fiB 
phase). Since the QCD phase transition is essentially the change of the degrees of freedom, 
we do not expect this to take place in continuum theory. In the strong copling lattice QCD, 
however, the baryon density is almost one in quark matter at low T, as discussed at T = 
in the previous section. In the maximum density case {pB = 1), all the sites are filled by Nc 
quarks and the entropy is zero. Thus the entropy density can be smaller in quark matter, 
while this is an artifact of finite spacing lattice. The same behavior is observed in staggered 



fermions 



11 



12], and it is found that the region with the positive slope boundary narrows 
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16|. 



or disappears with finite coupling effects in staggered fermion 

We note that in r 7^ 1 case, we have nonzero 714 condensate as shown in Fig. [6l where 
we show a and 7T4 as functions oi ps for several fixed T. This 714 condensate also undergoes 
a phase transition with the chiral transition. Since 714 is an auxiliary field for the operator 
iip'jAip, an absolute value of this vector condensate is deeply related to density. The point is 
that this condensate still remains for ps = 0- It is quite natural since the 0(1) renormalized 
imaginary chemical potential pi^ due to quantum effects can survive even if we set p^ in 
the minimal-doubling range to eliminate the 0(1/ a) imaginary chemical potential, as we 
discussed in Sec. IIIBi Thus, the physical imaginary chemical potential is given by the sum 
of the effective one /iim and a usual 0(1) parameter pi^ as Pi^'"^ = pim + Pim- It suggests 
that one possible way to control Pi^'^ in lattice QCD is to check the size of 714 condensate. 

The existence of pi^ also affects the ps dependence of critical lines: In Fig. [7] we change 
Ps and depict a three-dimensional chiral phase diagram for T, pb and p^ for r = 1. It 
shows that the critical line changes with p^ being varied. As far as we keep ps in the 
physical range, the 0(l/a) effective chemical potential for the physical flavors is cancelled. 
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FIG. 7: Three-dimensional chiral phase diagram for T, fis and fi^ for m = where /is runs within 
half of the physical range —3 < fi^ < ^f^ljl — 3. Green, red and purple lines show 2nd, 1st order 
transitions and tricritical point, respectively. The critical chemical potential at T = is explicitly 
given by (pT 



As discussed above, however, we still have 0(1) contribution as /iim- We can interpret that 
the dependence of the critical line on /i3 comes from the dependence of /iim on /i3 as /iim(/^3)- 
In this section we have obtained the finite- (T,/i) QCD phase diagram in the strong- 
coupling limit. To be precise, since the theory effectively contains the renormalized imagi- 
nary chemical potential as the KW artifact, it should be called the finite-(T',yU.Re)/^im ) QCD 
phase diagram. Anyhow, we have shown that we can apply the KW fermion to in-medium 
lattice QCD. We lastly discuss the real-type FCP fermions 
can also consider the real-type KW fermion, which loses 75 hermiticity. We can perform the 
strong-coupling QCD analysis for this type in a parallel way, and can derive the QCD phase 
diagram as long as the relevant parameter is set to the physical range. In the practical lattice 
QCD simulations, however, we should encounter a severe sign problem with the real-type 
FCP fermion even for zero-density cases. We need further study to judge its applicability 
to lattice QCD. 



26| . As shown in Sec. |TT1 we 
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IV. SUMMARY AND DISCUSSION 

In this paper we propose a new framework for investigating the two- flavor finite- (T,yu) 
QCD phase diagram. We show that the discrete symmetries of the Karsten-Wilczek (KW) 
fermion strongly suggest its applicability to the in-medium lattice QCD. To support our 
idea, we study the strong-coupling lattice QCD in the medium and derive the phase dia- 
gram of chiral symmetry for finite temperature and chemical potential. We have obtained 
the phase diagram with 1st, 2nd-order and crossover critical lines, which is qualitatively 
in agreement to results from the model study. We also argue that the additive chemical 
potential renormalization to the two flavors can be controlled by adjusting the parameter of 
the dimension-3 counterterm. 

In Sec. |TT] we review the Karsten-Wilczek fermion and its discrete symmetries. By the 
careful comparison between the KW fermion and naive lattice fermions with chemical po- 
tential, we find that both possess the same discrete symmetries, and that KW fermion is a 
proper formulation for finite- (T,/i) lattice QCD study. It is natural since the KW fermion 
decouples 14 doublers by assigning them 0(l/a) imaginary chemical potential, which also 
leads to the additive chemical potential renormalization for the rest 2 flavors in the interact- 
ing theory. We discuss that, in order to control this effective imaginary chemical potential, 
we need to keep a value of the relevant parameter fi^, within the minimal-doubling range 
in Fig. [21 It is more modest but similar tuning to the mass parameter tuning in Wilson 
fermion. One possible indicator of the parameter ranges is the pion spectrum: If /is is in the 
no-flavor range, there is no SSB of chiral symmetry and no massless pion. If fi^ gets into 
the six-flavor region, the number of pseudo Nambu-Goldstone bosons increases. In Sec. IIIII 
we perform the strong-coupling flnite-(T,yu) lattice QCD with the KW fermion. We derive 
the effective potential of a and tt^ as a function of temperature T, chemical potential /i, and 
/i3. As long as keeping fi^ within the physical range in Fig. [21 we successfully obtain the 
QCD chiral phase diagram, which is consistent with the phenomenological predictions: For 
high temperature or large chemical potential the chiral symmetry restores with chiral con- 
densate's disappearing. Our result strongly suggests that the Karsten-Wilczek fermion, or 
more generally flavored-chemical-potential fermions are useful for the 2-flavor lattice QCD 
with temperature and density. 

One potential problem for this formulation is that, even if we can set fi^ in the minimal- 
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doubling range, the rest 14 species with infinitely large chemical potential may contribute 
to the thermodynamical quantities. In such a case we need to subtract a divergent part 
properly. Further study is needed to figure out this problem. 

Finally we comment on the lattice QCD with this formulation in the small chemical 
potential limit. The minimal-doubling fermion breaks the spatial symmetry even for a zero 
chemical potential limit unless we fine-tune two more parameters to restore the hypercubic 



symmetry 



29(1 . We thus expect that this formulation could be less effective to describe the 



two-flavor QCD for smaller chemical potential although it works for the region near critical 
lines with sufficiently large chemical potential. 
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